function [lower_bounds, upper_bounds, best_wkde] = AcombinePI(bounds, test_ori, alpha)

    % 优化权重w
    SearchAgents_no=30;         % 种群数量，Number of search agents
    Max_iter=5;                % 最大迭代次数，用于计算a,a2。用于zeros函数生成一个1*10的矩阵
    dim=1;     
    lb_wkde = 0;
    ub_wkde = 1;
    
    fobj = @(w_kde) calc_conbinew(w_kde, alpha, test_ori, bounds);
    tic;
    [best_fun,best_wkde,cuve_f] = WOA(SearchAgents_no,Max_iter,lb_wkde,ub_wkde,dim,fobj);
    toc;
    
    w_btstrp = 1-best_wkde;
    best_lower_bounds = best_wkde*bounds(:,1);
    best_upper_bounds = best_wkde*bounds(:,2);
    best_lower_bounds = best_lower_bounds + w_btstrp*bounds(:,3);
    best_upper_bounds = best_upper_bounds + w_btstrp*bounds(:,4);
    lower_bounds = best_lower_bounds;
    upper_bounds = best_upper_bounds;
    lower_bounds(lower_bounds<0) = 0;
    upper_bounds(upper_bounds<0) = 0;
    
    % 计算最终的区间估计
    PICP = calc_PICP(test_ori, best_lower_bounds, best_upper_bounds);
    PINAW = calc_PINAW(test_ori, best_lower_bounds, best_upper_bounds);
    % 显示测试集的 PICP 和 PINAW 的结果
    fprintf('%d%%PICP: %.2f%%\n', (1-alpha)*100, PICP*100);
    fprintf('%d%%PINAW: %.4f\n', (1-alpha)*100, PINAW);

    % CWC = calc_CWC(PINAW,PICP,alpha,0.8);
    % fprintf('%d%%CWC: %.4f\n', (1-alpha)*100, CWC);

    AIS = calc_AIS(test_ori, best_lower_bounds, best_upper_bounds, alpha);
    fprintf('%d%%AIS: %.4f\n', (1-alpha)*100, AIS);
end